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Abstract 

A slave-boson functional-integral method has been developed for the Hubbard model 
with arbitrary, orbital degeneracy D. Its saddle-point mean- field theory is equivalent to 
the Gutzwiller approximation, as in the case of single-band Hubbard model. Our theory 
is apphed to the doubly degenerate {D = 2) model, whose paramagnetic state has been 
studied by numerical calculations. The effect of the exchange interaction on the metal- 
insulator (MI) transition is discussed. The critical interaction for the MI transition is 
analytically calculated as functions of orbital degeneracy and electron occupancy. 
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§1. Introduction 

The Hubbard model has been accepted as a model for a study of the strongly correlated 
systems such as transition metals, valence-mixing and high-Tc materials. The single-band 
Hubbard has been extensively investigated by using various methods. Gutzwiller^^ 
adopted a variational approach, employing the variational wave function. Since the exact 
evaluation of the ground-state energy is difficult, he used an additional approximation 
which is now called the Gutzwiller approximation (GA). A validity of the GA was in- 
vestigated by Monte-Carlo technique for finite-size clusters.^'^^ It is realized that the GA 
becomes a better approximation for the higher-dimensional system. ^'^^ The GA has been 
widely employed for various studies including magnetism and the metal-insulator (MI) 
transition. 

Kotliar and Ruckenstein (KR)^) proposed the slave-boson functional-integral method, 
whose saddle-point approximation is shown to be equivalent to the GA. The KR method 
has a wider applicability than the GA as follows: 

(1) within the saddle-point theory, we can deal with the complicated magnetic systems 
such as antiferromagnetic state besides paramagnetic and ferromagnetic states. 

(2) we can go beyond the saddle-point approximation, by including fluctuations around 
the saddle points. 

(3) we can extend the theory to finite temperatures. ^^'■^'''^ 

In contrast with the single-band Hubbard model, little studies based on the GA had 
been reported on the degenerate-band Hubbard model since its first application was made 
by Gutzwiller and Chao.^^^ In the last few years, however, several attempts-*^^"^^^ were 
reported of re-formulating and applying the GA to degenerate-band Hubbard model. Lu^^^ 
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obtained the analytical expression of the critical interaction for the MI transition. Okabe^'^ 
proposed a sophisticated method in calculating the band-narrowing factor, which is the 
most difficult part in applying the GA. Quite recently Biinemann and Weber^^^ discussed 
the first-order MI transition in the half-filled doubly degenerate Hubbard model when the 
exchange interaction is included. 

It is desirable to develop the slave-boson functional-integral theory for the degenerate- 
band Hubbard model, which is the purpose of the present paper. We have developed 
such a theory, adopting the slave-boson method of Dorin and Schlottman^^^ which was 
originally employed for the Anderson lattice model. As in the single-band Hubbard model, 
the saddle-point approximation to our slave-boson theory is equivalent to the GA in the 
degenerate-band Hubbard model, and our slave-boson theory has the above-mentioned 
advantages with the wider applicability than the GA. 

The paper is organized as follows: In the next §2, we present a slave-boson functional- 
integral formulation for the Hubbard model with the i^-fold orbital degeneracy. Our 
formalism is apphed to the doubly degenerate {D — 2) model, for which numerical cal- 
culations are performed and reported in §3. Section 4 is devoted to conclusion and sup- 
plementary discussion. In the Appendix we analytically calculate the critical interaction 
strength for the MI transition, in order to demonstrate the feasibility of our theory. 

§2. Formulation 

§2.1 Basic Equations 

We adopt the Hubbard model with D-iold orbital degeneracy whose Hamiltonian is 



given by 

-^"^^X^^^^ij ^ima^.im'a + 2 ^ E ^ mm' ^ima^imaC^im' a' ^irn' a' ■, (1) 
c mm' i (m,(T)7^(m',(T') 

where Cimc is an annihilation operator of an electron with an orbital index m and spin cr(=| 
, I) on the lattice site i and t^™' is the hopping integrals. As for the on-site interaction, 
^mm'i introduce three parameters: Uq for m — m', a ^ cr', Ui for m ^ m', a ^ o"', and 
U2 for m 7^' m', cr = cr'. We assume that the spin and orbital degrees of freedom are not 
coupled for simplicity. 

We employ the slave boson formulation of Dorin and Schlottman,^^^ introducing 2^^ 
bose annihilation operators: 

h^^^ (2) 

i;m\cri,m2(T2,....,me,cre,i \ / 

which project to the configuration of I electrons with pairs of orbital and spin indices 
{ma}. In eq. (2) all the index pairs are different because of the Pauli principles. 
We introduce the following product operations:^^^ 

(1) A full contraction is defined as 

y-'i '-'i J / J ^i;mi(Ti,m2(T2,----,mfa( i;mi<Ti,m2(T2,----,m((T(' \ J 

miai ,m2cr2,....,m£a£ 

(2) A partial contraction with a fixed pair of (m„cr„) for £ > 1 is defined as 

\ I I Ji^nO'n — / J i;miai,m2cr2,....,m(a( i;miai,m2cr2,....,meai' 

miai,m2<72,----,m„-i<Jn-i,mn+icrm+ii,----,me<7e 

(4) 

With the use of these operations the completeness of the slave boson operators is given 

by 

2D 

E(bf ^ • ) = 1, (5) 



while the correspondence between bosons and fermions is expressed by 

2D ^^^^ 

The partition function of the model given by eq. (1) is expressed as a functional 
integral over coherent states of Fermi and Bose fields. The constraints given by eqs. (5) 
and (6) are incorporated with Lagrange multipliers, A^^^ and A^^^^ to get 

Z^l Uj^,Db('^ J m(^) / DX(^^ exp[- £ dr {Lf{r) + L,(r))], (7) 

where the actions, Lf and Lf„ are given by 

Lf (-r) = E E [4ma {d/dr + ASi) Q^^] + E E E ""LacLaCjm'aZjma, (8) 

i ma c ij mm' 



i=0 mcr e=i 



+ *o, (9) 



2D 

^ = "V "V \" \^ TTTa' M] M /-|A\ 

/ J / J / J / J mm' i;miai,m2cr2,....,m(a( i;miai,m2(T2,----,m(a('> \ J 

i 1=2 micn,m2<T2,-;meat {ma,m'a') 

the fourth summation in eq. (10) being made over a pair of indices {ma^m'a') with 
(mcr) 7^ {m'a') in the configuration: {micri, m2cr2, m^cr^} occupied by I electrons. It 
is noted that when C/q = C/i = C/2 = C/, we get $0 = Ei 2 ^"2 (bf^^ • b?) where 
Cj^ — -^l/A;! (£ — /c)!. Although a choice of the factor, Zim^, is not unique, the expression 
given by 

2D 

Z,ma-{l-n,mar/'T.{bt'^'-^^) (H) 

~~ \ ' ma 

l=\ 

is shown to yield the correct result in the limit of vanishing interactions.^^-* 

We employ the static approximation for boson fields. Furthermore we adopt the 
following change of variables: 

^im = (Aj-J| + A-^^)/2, i,im = -(A-^-f - AjJj/2, (12) 



^im — < T^iml > + < Uim^ >, TTlim — < ''T-imt > ~ < '^imt (13) 

where the brackets, < >, stand for the average. Thus the expression of the partition 
function becomes 

Z = J U'j22Db''^^ j Dv j Di j Dn j Dm exp (14) 



where 



e ^'^ — exp ( —(3 



i m 



Tr exp{-pHeff). (15) 



The effective one-electron Hamiltonian, i^e//, is given by 



Heff — Yl ^ima^jm'atij '^ima'^jrn' a + '^^'^^'^^iyi'in (^Cim) C^ma^imaj (1^) 

c ij mm' aim 



z- -2(2-n- r^/^n"^/^ 
where and Pi^a ^-re given by 



2D 

Pima + VPi^ + XI (^f ■ 



(17) 



2D 

= (bf . bf ) = 1 - E^- + EE [(^ - 1)/^] (bP • (18) 

m ma (,=2 

2D 

Pima = (bf ^ • ^)^(, = (rij^ + amim)/^ - E(bf ^ • hf^)ma- (19) 

It should be noted that the functional integral given by eqs. (14)- (19) is performed over 
variables of ^, i/, n, m, and h^^^ with £ > 2. Our expression has physically more transparent 
meaning than an original static approximation in which the functional integral is carried 
out over the variables of A^^\ A^^-' and h'^^'' with £ > 0, although both the approaches 
are equivalent. The expression for the functional integral given by eqs. (14)-(19) is a 
simple generalization of the single-band model to the degenerated-band model. When 



7 

all the integration variables are replaced by their saddle-point values, we get the mean- 
field approximation, which is expected to be equivalent with the GA theory, ^^"^^^ related 
discussions being given in §4. 



§2.2 Doubly degenerate band 

We apply our theory developed so far to the doubly degenerate Hubbard model. After 
a simple calculation, the expressions for the functional integral given by eqs. (14)-(19) 
become 

Z = J DdoJ DdiJ Dd^J Dt^jDf J DC j Dm j Dv j Dn exp[-/3($o + + ^2)], 

(20) 

$0 = E Wodi^ + 2C/ida + U2{di^ + d^^) + 2(C/o + C/i + U2)m + + h)] , (21) 

i 

^1 = 1] [iimmim + ^im)nim\ , (22) 

im 

$2 = / def{e){-l/T:)Im Tr lnG{e). (23) 

where f{e) is the Fermi distribution function and n Fermi level. The one-particle Green 
function, G(£), is expressed in terms of the effective Hamitonian, i^e//, for the system 
under the charge {yim) and exchange (^i^) fields, as 

G{e) = (£ - H,if)-\ (24) 

with 



-^eff ^ima^ij ^jm'a '^imcr^jm'cr + ^^^^^^(^im '^dm) ^imu'-'ima i (^5) 

f ij mm' aim 



^JPiaWei + ^/di^) + {y/dio + Vdii){VPi-a + Vk^ + ^/t~{y/d^^^ + ^1 fi. 



(n-im + (ynriimY/'^ (2 - - (7mi^)V2 



,(26) 



Ci = 1 - 2nim + 2dio + 2dii + di^ + + A{ti^ + t^^) + 3/^, (27) 

Pia = (ni,„ + amim)/2 - {dio + dn + dia) - {2tia + U-a) - fi, (28) 

where tj^- and stand for the i = 3 and 4 components of (bf'' • bf^)^^., respectively. As 
for £ — 2,we take into account the three kinds of configurations: dio for a pair of electrons 
on the same orbital with opposite spin, dn on the different orbital with opposite spin, and 
dia- on the different orbital with same spin a. 

The exphcit form of the Green function, G{e), depends on the electronic and magnetic 
structures of a system to be investigated. Since the effective transfer integrals in eq. (25) 
is expressed as a product form: Zimatij^' Zjm'a, we can express the one-electron Green 
function in terms of the locators defined by^^-* 

where Vima — {^imaT, with which the band-narrowing factor, q^^' , is expressed as 

The mean-field free energy can be obtained from the saddle-point values of the in- 
tegration variables for which the variational conditions yield the following simultaneous 
equations: 

cr 

fTlim — ^ ^ < Cj^^^Ciraa ^i (32) 
- i^im + X! ^irrxa {9^^^/ dfli^) = 0, (33) 
6m + Rima{drima/dmi^) = 0, (34) 



2t^0 + ^Rima{dri,na/ddiQ) = 0, 
ma 

2t/l + ^Rima{dri,na/ddii) = 0, 
mcr 

mcr' 

2{Uo + Ui + U2) + J2 Rima'{dn^a,/dti„) = 0, 
ma' 

2{Uo + Ui + U2)+Yl RimaidUma/dfi) = 0, 
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(35) 
(36) 
(37) 
(38) 
(39) 



where Ri^a = d^2/drima- 

In the remainder of this paper we consider only the paramagnetic states at T = K, 
for which we get 

$2 = pde {e-fx) {-l/n)Im TrG{e), 



q = r — 



{e-v - qsk) ' 
£o(?T') = ^Ra, 

a 

( V'^ + V'^ + y/di) (y/p + y/t) + y/ep+ y/tff ^ 



(40) 
(41) 
(42) 



with 



n(2 - n) 



e = 1 - 2n + 2{do + di + da) + St + 3/, 



p^n/2- {do + di + da) -^t-f. 



(43) 



(44) 



(45) 



Here d2 — d^, t — and R^ — Rima al, subscripts i and m being omitted; Ek is 
the Fourier transform of the transfer integrals and £o(^) is the ground-state energy per 
sub-band for C/q = C^i = C/2 = as a function of n, the number of electrons per sub-band. 
The ground-state energy per site is then given by 



E = 2qeo{n) + 2Uodo + 2Uidi + 2U2d2 + 2{Uo + Ui + U2){2t + /), 



(46) 
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which is identical with that obtained by Bliimenann and Weber. ^^-^ Employing the 
Gutzwiller-type wave function, they obtained the ground-state energy given by eq. (46), 
from which the optimum values of the occupancies are determined by the variational 
conditions for E given by 

Uo + eo(dq/ddo)^0, (47) 
Ui + eo{dq/ddi) = 0, (48) 
U2 + eo{dq/dd2)^0, (49) 
2{Uo + Ui + U2) + eo{dq/dt) = 0, (50) 
{Uo + t/i + U2) + eo{dq/df) = 0. (51) 
It is easy to see that eqs. (35)-(39) reduce to eqs. (47)-(51) in the paramagnetic states. 

§3. Numerical Calculations 

Numerical calculations have been performed for D = 2 model with the quarter-filled 
(n = 0.5) and half-filled {n — 1) bands, which are considered to be most interesting. 

We firstly show the calculated results by taking Uq — Ui — U2 — U, for which do, di and 
d2 are equivalent: do = di = d2 = d. The ?7-dependence of the occupancies for n — 0.5 
is shown in Fig. 1. At C/ = we get e = 0.3164, p = 0.1055, d = 0.0352, i = 0.0117 
and / = 0.0039. When U value is increased, e, t and / decrease while p increases. In 
particular, t and / become at least ten times smaller than e, p or d ai U / \ eo{n) \ > 3. The 
q factor also decreases and vanishes above [4 which stands for the critical interaction for 
the metal-insulator (MI) transition. Near the MI transition, q behaves as g oc {Uc — U) 
as in the single-band model. We get Uc — 13.2 | £o(^) I, which agrees with the previous 
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results. ^^'^^^ Above Uc only p remains finite with a value of 0.25. The U dependence of the 
ground-state energy is shown in Fig. 2, where the energy calculated by the Hartree-Fock 
approximation (HFA): 

En^ = 2£o(n) + (1/2) (C/q + C/i + C/2) n\ (52) 

is also plotted for a comparison. The ground-state energy in the GA is always lower than 
-Ehf as expected. 

Figure 3 shows the U dependence of the occupancies for half-filled (n = 1.0) case, in 
which the relation: 

e = /, P^t, (53) 

holds by the electron-hole symmetry. When U value is increased, only d increases while 
other occupancies decrease. The band-narrowing factor, g, vanishes above [4, which is 
12.0 I £o(^) I in agreement with the result of refs. 19 and 21. The double occupancy d 
persists even above Uc with a value of 0.167, which is in contrast with the quarter-filled 
case where only single occupancy p survives. We might expect that in the limit of infinite 
[/, any states with more than doubly occupied state vanish. This seems not the case 
for the half-filled D — 2 degenerate band, where doubly occupied state may survive by 
making a compromise with the constraints given by the completeness of bosons (eq. (5)) 
and the electron-hole symmetry (eq. (53)). This situation is changed when the exchange 
interaction is introduced, as will be discussed shortly. Fig. 4 shows the ground-state 
energy of the half-filled case as a function of U . Again the energy obtained in our theory 
is lower than that in the HFA. 

Next we take into account the exchange interaction, J, with which our three interaction 
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parameters are expressed as 

Uo^U, Ui^U' - J -2J, U2^U' - J -3J, (54) 

where we employ the relation: U — U' — 2J required for the rotational symmetry, U 
and U' being the intra- and inter-band Coulomb interactions, respectively. 

Figures 5 and 6 show the U dependence of the occupancies and the ground-state energy, 
respectively, for the quarter- filled band with J/U = 0.1. The general behavior shown in 
Figs. 5 and 6 is similar to that in Figs. 1 and 2, respectively, except the fact that the 
degeneracy among doubly occupied states is removed with an inclusion of the exchange 
interaction, J. Because Uq > Ui > U2 for finite J, we get do < di < ^2- We realize that 
an inclusion of the exchange interaction also enhances the value of Uc- Ud \ e{n) \ = 13.2, 
16.4, and 24.8 for J/U — 0.0, 0.1 and 0.2, respectively. The MI transition is of the second 
order in all the quarter-filled cases investigated. 

On the contrary, the effect of exchange interaction is quite different in the case of 
half-filled band. The MI transition becomes the first-order transition as is shown in Fig. 
7, where the U dependence of the occupancies with J/U — 0.1 for n — 1.0 is plotted. 
This first-order MI transition has been recently pointed out by Biinemann an Weber. A 
change from the second-order to first-order MI transition is induced with an infinitesimally 
small inclusion of J. When J = 0, all doubly occupies states are equivalent and remain 
finite even at t/ > U^ as was shown in Fig. 3. For a finite J, however, d2 is the 
most favorable state among the three doubly occupied states, and only d2 survives above 
U — 6.0, where neither do nor di remain finite. This first-order transition is realized in 
the U dependence of the ground-state energy shown in Fig. 8, where the E — U curve 
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of the metallic state does not tangentially intersect that of the insulating state at the MI 
transition, in contrast with the case of the second-order transition shown in Fig. 2, 4 and 
6. 

§4. Conclusion and Discussion 

To summarize, we have developed the slave-boson functional-integral theory for the 
Hubbard model with the orbital degeneracy D, employing the slave-boson method of 
Dorin and Schllotman^^'' who originally applied it to the Anderson model. Our slave- 
boson mean-field theory has been used for a study of the doubly degenerate model, for 
which numerical calculations have been performed. We have shown that our slave-boson 
mean-field theory is equivalent with the GA for D = 2.^^^ It is shown in the Appendix 
that our mean-field theory and Lu's^^^ yield the identical critical interaction of the MI 
transition for arbitrary D. This suggests that the equivalence between the slave-boson 
mean-field theory and the GA^^"^^) generally holds for any degeneracy D, although it has 
not been rigorously proved yet. 

We expect that our slave-boson mean-field theory is a convenient and useful method in 
studying strongly correlated systems described by the degenerate-band Hubbard model. 
With the use of the Green function method, we can deal with the complicated magnetic 
states within the mean-field approximation. It would be indispensable to investigate the 
antiferromagnetic states^^^ for a deeper understanding of the MI transition of the D — 2 
Hubbard model. Our slave-boson mean-field theory enables us to investigate the occupa- 
tion versus interaction phase diagram of the degenerate Hubbard model, which was previ- 
ously studied within the Hartree-Fock approximation. Furthermore, we can take account 
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of fluctuations around the saddle points, going beyond the mean-field approximation. It is 
interesting to extend our theory such as to include the effect of spin fluctuations at finite 
temperatures in order to discuss the finite-temperature magnetism of transition metals. 
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§ Appendix 

As one of applications of our formalism, we calculate the critical interaction, Uc, for 
the MI transition of the Hubbard with the orbital degeneracy, D, in the paramagnetic 
states. We assume that near the second-order MI transition point, fluctuations of the 
electron occupancy around its average N (integer) are small. Then we take into only 
three possible states: N — 1,N and + 1 states, whose occupancies are expressed by e,p 
and d, respectively, after Lu.-*^^^ 

The constraints given by eqs. (5) and (6) become 



y2D 



2D 



(A-l) 



CJTS,' e + Cl^S,' p + C-jT-' d = n/2, 



2D-1 



(A-2) 



where n (= N/D) is the number of electrons per sub-band and Cf [= i\/k\ {i — k)\ ] is 
zero when either £, k, or {£ — k) is negative as convention. The ground-state energy when 
Uq = Ui = U2 = U m. eq. (1) is given by 



E^qD £o(n) + U {cl^^^.C^' e + Cl^'C^ p + Cl^i.C^^' d 



(A-3) 



where e and p are given as a function of d as 



q = [4/n(2 - n)] + ^^^^^ Vp A 



(A-4) 



e = (Ai/A) d, 



(A-5) 



p = (A2/A) - (A3/A) d. 
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with 



N-2 1 



(A-7) 



A _ r<2D-1^2D ^2£)-1^2D 



(A-8) 



A, = (n/2) Cfli, - Ci,%\ 



(A-9) 



A — /^2Z)-1^2D r^2D-lr^2D 



(A-10) 



Taking the variational condition for E with respect to d, we obtain the critical interaction, 



Uc, given by 



D 



£o(n) I N{2D - N) 



^N{2D - AT + 1) + ^{2D -N){N+ 1) 



(A-11) 



which is identical with the result obtained by Lu.^^^ Equation (A-11) is symmetric with 
respect to the interchange of the variable: N ^ {2D — N) as it should be. In the limit of 
infinite D and finite N, eq. (A-11) becomes 



lim 

D-*oo 



Ur 



D I £o("') 



N 



(A-12) 



Equation (A-12) is in agreement with the result of ref. 26, in which the fermion problem 
is treated as the boson problem using the GA. Some calculations using Eq. (A-11) are 
shown in Fig. 9. 
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Figure Captions 

Fig. 1 The occupancies, e,p, d, t and /, and the band narrowing factor, 5, as a function of 
U for the quarter-filled {n — 0.5) case with J = 0, the results of t and / being multiplied 
by a factor of ten. 

Fig. 2 The ground-state energy for n — 0.5 and J = in the GA (sohd curve) and the 
HFA (dashed curve), the arrow denoting the MI transition. 

Fig. 3 The occupancies, e,p, d, t and /, and the band narrowing factor, g, as a function 
of U for the half-filled (n = 1.0) case with J = 0. 

Fig. 4 The ground-state energy for n = 1.0 and J = in the GA (sohd curve) and the 
HFA (dashed curve), the arrow denoting the MI transition. 

Fig. 5 The occupancies and the band narrowing factor, q, as a function of U for the 
quarter- filled {n — 1/2) case with J/U = 0.1, the results of t and / being multiplied by a 
factor of ten. 

Fig. 6 The ground-state energy for n — 0.5 and J/C/ = 0.1 in the GA (solid curve) and 
the HFA (dashed curve), the arrow denoting the Ml transition. 

Fig. 7 The occupancies and the band narrowing factor, q', as a function of U for the 
half-filled (n = 1.0) case with J/U = 0.1. 
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Fig. 8 The ground-state energy for n = 1.0 with J/U = 0.1 in the GA (solid curve) 
and the HFA (dashed curve) , the dotted curve showing the energy of the insulating state 
{q — 0) and the arrow denoting the MI transition. 

Fig. 9 The critical interaction for the MI transition, Uc/ \ So{ri) |, against the electron 
occupancy per sub-band {n — N/ D) of the Hubbard model with electrons and D-fold 
orbital degeneracy. Only points shown by filled circles are meaningful, dashed curves 
being drawn using eq. (A- 11). 



